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Abstract 

Standard results for the fluctuations of thermodynamic quantities are derived under the assump- 
tion of sampling identical systems that are in different, not fully equilibrated states. These results 
apply to fluctuations with time in a particular macroscopic body and can be traced to the fluc- 
tuations of thermal occupancy. When many identically prepared - but not identical - systems are 
studied, mesoscopic fluctuations due to variations from sample to sample contribute to the fluctu- 
ations of thermodynamic quantities. We study the combined effect of mesoscopic fluctuations and 
fluctuations of thermal occupancy. In particular, we evaluate the total particle number and specific 
heat fluctuations in a two-dimensional, non-interacting electron gas in classically integrable and 
chaotic circumstances. 
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I. INTRODUCTION 



The state of thermal equihbrium of a closed system is achieved when its entropy is at 
maximum. In this state, the thermodynamic quantities characterizing macroscopic bodies 
that comprise it are at their mean values (thermodynamic averages) . [l| Thermodynamic 
fluctuations are the deviations of these quantities from their mean values and are described 
by the Einstein formula for the probability density [2I: 

w {x) = const exp [S (x)] 

where S (x) is the entropy and x = {xi,X2, ■ ■ ■ ,Xn} are the thermodynamic quantities. 
Ultimately, these fluctuations can be related to the fluctuations of thermal occupancy. For 
instance, thermal occupancy of quantum gases can be derived starting with a generally non- 
equilibrium state, followed by maximization of entropy subject to constraints (conserved 
quantities) . 

This approach to fluctuations, however, neglects sample-specific variations: it assumes 
that sampling occurs in the ensemble of identical systems that are in different states, slightly 
away from the equilibrium.!^ On the other hand, it has become abundantly clear over the 
last half-century that sample-specific variations are an important source of fluctuations, as 
witnessed by developments in mesoscopic physics. These must be included when thermody- 
namic properties of many identically prepared, yet not identical, systems are studied. In this 
work we investigate the combined contribution of mesoscopic fluctuations and fluctuations 
of thermal occupancy to fluctuations of thermodynamic quantities. 

We use a simple formalism that gives fluctuations of thermodynamic quantities as a sum 
of two terms: the first due to thermal occupancy fluctuations and the second due to sample- 
specific fluctuations. The latter is expressed in terms of the correlation function of the 
density of energy levels. The results are illustrated for the fluctuations of the number of 
particles and of the specific heat in systems with the fixed chemical potential and volume. 
As a model system, we use the degenerate, two-dimensional, non-interacting electron gas 
with level statics corresponding to classically chaotic and integrable motions respectively. 



2 



II. FORMALISM 



In the most general form, the value of a thermodynamic quantity, whose density is ex- 
pressed as a function of energy, can be computed as the integral over the spectrum 

G = I dep{e)f{e)g{e) (1) 

where p (e) level density and / (e) is the thermal occupancy. Under an obvious assumption 
that p and / are uncorrelated variables, we find the following expressions for the mean and 
the variance 

G = [ dep{e)J{e)g{e) (2) 



5G^ = 5GI + 5G) (3) 




6Gl= deide26p{ei)6p{e2)f{ei)f{e2)g{ei)g{e2) (4) 



SGj = j j deide2p{ei)p{e2)Sf{ei)6f{e2)g{ei)g{e2) (5) 
where the over-bar stands for a mean value. Using 

where p is the chemical potential and T is temperature, we find 

(f « 

Combining eqs. ([2]), ([5]) and (jTj), we find, as expected from SGj = [dG/dN) SN + 
{OG/dT) 6T, 



dp J \dp J ^ \dT 
where G is the mean specific heat and we used the well-known results from the theory of 

a 

thermodynamic fiuctuations[2j 

dp' C 



SN] = T— , 5T2 = ■= (10) 



to obtain the second equality. Applying this result to the fluctuations of the specific heat, 
for instance, we find 
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III. FLUCTUATIONS IN A 2D ELECTRON GAS 



We apply the above results to a non-interacting, degenerate, two-dimensional electron 
gas occupying area A. In this case 

, , ON mA ^1 

is a constant inverse average level spacing and m is the electron mass. Using standard low- 
temperature (Sommerfeld) expansion, we find ^-Sp = — Texp (— ei^/T), where ep is the 
Fermi energy, and 
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SN] = -, 6Cj = C^ -- (13) 



Next, we proceed to evaluate SN^ and SC^ for the classically chaotic and integrable circum- 



stances respectively and compare those to SN"^ and 6Cj. 

To simplify calculation, we assume a unitary ensemble (broken time reversal symmetry) 
in the classically chaotic case ^ 



6p (e) 6p {e + uj) = 1 (14) 



6{uj) sin2(7rcj/A) 
A nu 

We also use a simplified ansatz in the classically integrable casejsl 



. . ^ . . , — V '^(^) sin(27r^/E^) 

6p{e)6p{e + u) = — 15 

Here is the the energy scale corresponding to the shortest classical periodic orbit. In 
a square well, for instance, Em = VnFK. The simplified ansatz neglects large oscillations 
of the density correlation function for e » E^,-^ These should be of little consequence 
at finite temperatures since harmonics with incommensurate frequencies will be mixed it to 
wash out the effect of such oscillations. Still, unlike in the chaotic case, the correlation 
function depends on e via E^- However, due to relevant e's being restricted to the interval 
of order T around ei?, we can replace e —>■ ep 

Em ~ ^/tvefA (16) 

With this substitution, it will be shown below that for T ^ Em the result for integrable 

for T > Em, the full oscillatory 
should be of little consequence 



systems is due to the delta function term in ( 1151) . Converse 
behavior of the level correlation function discussed in Ref. 



since their scale is defined by Em- This reasoning is ultimately confirmed by numerical 



evaluation using the full form of the correlation functionpl] (with substitution e ep), 



which yields results that are extremely close to the analytical results for ^A''^ obtained 
below using (ITSll . 

To eliminate temperature-independent divergencies, we evaluate 



''"^ - f (17) 



where we used that / de2Sp{ei) Sp{e2) = 0. Changing variables to e = (ei + £2) /2, a; 
£2 — £1 and integrating on e, we find 



dSN^ 
dT 

where 



/OO 
duj5p (e) 6p (e + uj)H (^) (18) 



H(x) = 



X {X 

— csch — 
2 V2 



1 2 



(19) 



We first substitute dHj), ( |T5l) and (fT9l) into (|T71) and evaluate the integrals analytically. 
Further integration on T gives 

and 



5iV2 = 5Ar2 - + coth -— - (21) 

for classically chaotic and integrable systems respectively, with given by f|T6l) and zero- 
temperature fluctuations given by 

imf) = r r de^de25p{e^)5p{e2) (22) 

V /T=0 Jq Jq 

The latter is simply the variance of the number of levels over the Fermi sea and we use the 
results for the level number fluctuations from Ref. 5|] to find 



E.. 



and 

in,T=0 A 

for classically chaotic and integrable systems respectively 




SV| , , (24) 
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FIG. 1: Function F (x) given by eq. (j3ip 



Combining now with f|T3l) . we find for tlie total particle number fluctuation SN"^ = 6N'^ + 



V Jch V PJch,T=0 27r2 27T^T/A J ^ ' 

(26) 



and 



T/A - log (47r2T/A) /27r2, T > A/27r2 

m,T=0 A \ Era ) 27r2A 

27r2TV3E„A, T < ^„/27r2 
T/A - E^/27r2A, r>E„/27r2 



5iV^).^ - (5iV.^)__. = x-^h ( ^ ) - (27) 

(28) 



for classically chaotic and integrable systems respectively. Notice, that for temperatures 
below the energy scale where level rigidity sets injj],j5| - A for classically chaotic and for 
classically integrable systems - the main effect of temperature is to reduce sample specific 
fluctuations of the number of particles relative to the zero-temperature limit. Furthermore, 
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FIG. 2: ((^Cp) vs. T/A, with limiting cases given by eq. ()32p shown as thin straight hnes. 



the temperature-dependent part of such fluctuations largely cancels the fluctuations ( fT3|) due 
to thermal occupancy. Incidentally, since in this limit the temperature-dependent sample- 
specific fluctuations are dominated by the 5- function term in ( IT^ and (fT5i) . the temperature- 
dependent particle number fluctuation can be found directly from 
d r°° , , -^^^^^^ , Id 



d5Nl 



deide2Sp (ei) Sp (£2)/ (^i) / (^2) 



1 

'a 



(29) 



dT dT Jo Jo - r V r V ... V X, . V AdTjQ 

confirming the results in fl2Ul) and fl^Tl) . Conversely, for temperatures above the level rigidity 
scale, fluctuations due to thermal occupancy dominate. 

Fluctuations of the specific heat (heat capacity) can be found fromQ] 

-dj (er) dj (62) 



= J J delde26p{el)6p{e2)- 



^T dT 



where 



Fix) = 



120 



X 



{Utt^ + x') coth j - 4 (Tvr^ + Stt'x') 



tJ J du;Sp{e)6p{e + u;)F(^^) (30) 

(31) 



csch 



and is plotted in Fig. ([T]). As a result we find the limiting behavior of the fluctuations 

aT/A, r<A 
'^^ ^ b{2n^y\T>A 



(32) 



0.030 - 



0.025 



0.020 



0.015 



0.010 - 



0.005 -/ 



0.5 



1.0 



1.5 



2.0 



T 



FIG. 3: ((5Cp)^^ vs. T/Em, with limiting behavior (j34p shown as thin straight Une. 



where the constants a and b are given by 
a = F (0) 0.997, b = 

and 



dx 



[F (0) - F (x)] 



1.067 



x^ 



(33a) 



(34) 



for classically chaotic and integrable systems respectively. Clearly, the specific heat due to 
sample-to-sample variations in integrable systems has, over the large range of temperatures, 
the same functional dependence as that due to the fiuctuations of thermal occupancy. 

We show the plots of [SCp) and [SCp).^ vs. T in Figs. ([2]) and ([3]) respectively. Some 
of the curvature in these plots can be traced to that in Fig. ([1]). Small bumps and dips 
on the line in Fig. ([3]) is due to poor numerical convergence of rapidly oscillating integrals. 
For T > Em, [5Cp).^ shows faster than exponential decay, as expected from integration of 
an exponentially decaying function (l3Ti) with a function rapidly oscillating around zero. In 
contrast with chaotic systems, in integrable systems the decay of fiuctuations of the specific 
heat at high temperature is consistent with the greater level rigidity at larger energy scales. 
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as discussed in Ref.la 



IV. CONCLUSIONS 

We showed that the fluctuations of thermodynamic quantities can be separated into 
those that are due to sample-to-sample (mesoscopic) fluctuations and those due to thermal 
occupancy fluctuations. For a 2D non-interacting degenerate electron gas, both effects were 
included in evaluation of the particle number and specific heat fluctuations for classically 
chaotic and integrable cases. In the integrable case, the specific heat fluctuations due to 
either effect are comparable over a wide range of temperatures. For the particle number 
fluctuation, the temperature-dependent part of mesoscopic fluctuations is negative and, for 
temperatures below the energy rigidity scale, largely cancels the fluctuation due to thermal 
occupancy fluctuations. For higher temperatures, the latter dominates the temperature- 
dependent part of the fluctuation. 
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